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Abstract: We consider tilings of quadriculated regions by dominoes and of tri- 
angulated regions by lozenges. We present an overview of results concerning 
tileability, enumeration and the structure of the space of tilings. 



Introduction 

Let R be a finite juxtaposition of unit squares in the plane. We consider tilings of R by 
dominoes, 2x1 rectangles. Alternatively, if R is obtained from unit equilateral triangles, 
we try to tile R by lozenges (or calissons), unions of two adjacent triangles. The following 
questions are natural: 

• Can R be tiled? 

• How many tilings exist? 

A next set of questions is suggested by the concept of Hip. We perform a flip by lifting 
two dominoes forming a square or three lozenges forming a hexagon and placing them back 
after a rotation of 90° or 60°. 

• Do all tilings admit flips? 

• Are any two tilings joined by a sequence of flips? 

• What is the distance (the minimal number of flips) between two tilings? 

• Which flips decrease the distance to a given tiling? 

• Are there closed geodesies of flips in the space of tilings? 

Furthermore, as we shall see, there is a natural bipartition on the space of domino 
tilings so that flips connect tilings in opposite classes. 

• Is the number of tilings in each class the same? 

For simply connected regions, these questions are easier. At the other end of the 
spectrum, they still make sense for quadriculated or triangulated surfaces. 

Section 1 deals with the problem of tileability; the main references are [CL] and [T] . In 
section 2, we discuss some techniques for counting tilings, from the beginning of the century 
with [M] to recent contributions by mathematical physicists ([K] and [LL]). Finally, section 
3 deals with flips and related issues; the results were originally presented in [STCR]. 
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1. Existence of tilings 



Given a finite set of types of tiles, can we use them to tile the plane? Rather surpris- 
ingly, this is an undecidable question ( [B] ) . The corresponding problem for finite regions is 
NP-complete ([GJ]). Colouring arguments are among the simplest to provide necessary con- 
ditions for tileability. Both the triangulated and the quadriculated plane are bicolourable 
and a region must have the same number of triangles (or squares) of each colour to be 
tiled by lozenges (or dominoes). Conway and Lagarias ([CL]) discovered a method to 
obtain strong necessary conditions for tileability. Following Thurston ([T]), we present 
their method for simply connected quadriculated or triangulated regions, the tiles being 
dominoes or lozenges. 
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Figure 1.1 



We start with lozenges. In the triangulated plane, we may label the oriented edges a, 
b and c as in figure 1.1. a and paths along edges yield words in the letters a, b and c. Notice 
that a path is closed iff the total numbers of a's, 6's and c's are equal, i.e., iff its word equals 
the identity in G± = Z 2 , given in terms of generators and relations as G\ = (a,b,c\abc = 
acb = e). Consider the group G = Z 3 = (a,b, c\[a, b] = [a,c] = [b,c] = e) (as usual, 
[g, h] = g~ 1 h~ 1 gh): the relations are the boundaries of the three possible orientations of 
lozenges (i.e., position up to translation). The boundary of a simply connected region R 
yields a word w and thus an element g G G: a tiling of G gives a way to write w as a 
juxtaposition of conjugates of boundaries of tiles (see figure 1.2 or [CL] for an actual proof) 
and thus g = e. 




Figure 1.2 



These group-theoretical considerations produce Conway-Lagarias's condition for tile- 
ability of simply connected regions. We state their result for triangulated (or quadricu- 
lated) regions: 
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Theorem 1.1: Let R be a triangulated (or quadriculated) simply connected bounded 
region and E a finite set of (connected and simply connected) tile shapes. Let G be the 
group with generators corresponding to edges and relations given by boundaries of elements 
of E . Then a necessary condition for the tileability of R by translates of elements of E is 
that the word induced by the boundary of R is trivial in G. 

For lozenges, this condition reduces to the simple colour condition. More interesting 
applications of the group theoretic conditions yield the following results ([CL]). Let T n be 
the triangular region of size n in the hexagonal lattice (fig 1.3. a) and L 3 be the tribone, 
the juxtaposition of three hexagons in a line (fig 1.3.b). 

Theorem 1.2: The triangular region T n can be tiled by copies ofT 2 iff n = 0, 2, 9 or 11 
(mod 12). Also, T n cannot be tiled by congruent copies of L 3 . 




(a) (b) 

Figure 1.3 



David and Tomei noticed that a lozenge tiling (as in figure 1.4. a) may literally be 
seen as a picture of a bunch of cubes stored in a big box (as in figure 1.4.b). Align the 
sides of the box with the coordinate axes so that vertices of triangles receive integer x,y,z 
coordinates. In figure 1.4, for example, the box has size 5x5x5 and there are 66 cubes 
in it. The figure is an orthogonal projection of the three dimensional configuration on the 
plane x + y + z = 0. If we think of gravity as pointing "downwards" to (—1, —1, —1), the 
cubes form a gravitationally stable pile iff they represent a tiling. By projecting on the 
coordinate planes, we immediately obtain for hexagonal regions the following invariance 
result ([DaT]): 

Theorem 1.3: Let R be an arbitrary triangulated planar region. Then the numbers 
n xy, n xz , n yz of lozenges in each orientation do not depend on the tiling. 

The visual argument works for simply connected regions; for general regions we might 
use height sections (to be introduced later) but we present an elementary proof which 
appears to be folklore. Consider vectors going from the origin of the plane to the centre 
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Figure 1.4 



of each white triangle and from the centre of each black triangle to the origin and add 
them up obtaining a vector v. Let w xy be the vector joining the centres of triangles in 
a xy-lozenge, from black to white; define w xz and w yz similarly. Given a tiling, we may 
group pairs of vectors associated to the same lozenge to obtain: 

Since n xy + n xz + n yz does not depend on the tiling and w xy — w xz and w xy — w yz are 
linearly independent, the theorem follows. 

Thurston obtained the same three dimensional interpretation as a spinoff of the group- 
theoretical construction. Consider the Cayley diagram V for G and its natural projection 
onto Ti, the Cayley diagram of its quotient Gi (Ti is, of course, the triangulated plane): 
a tiling of a simply connected region lives in Ti but may be lifted to V (since each tile 
can). Imbed the Cayley graph T of G = Z 3 in R 3 (in the obvious way) and Ti in the plane 
x + y + z = 0: the projection of V onto Ti is the restriction of the orthogonal projection 
onto this plane. This lifting assigns to each vertex of a triangle the same coordinates as 
the more visual procedure represented in figure 1.4. 

The extra information obtained by lifting is x + y + z, the height with respect to 
the plane x + y + z = 0: we call the function which assigns to each vertex the sum of 
its coordinates a height function, as shown in figure 1.5. a. Of course, there are other 
expressions which convey the 3D information equally well but, for height functions of 
lozenge tilings, we always take x + y + z. Height functions are therefore well defined 
for a tiling by lozenges of a simply connected region up to an additive constant. For a 
given tiling, a simple local rule obtains the height function: when following an edge at 
the border of a tile, add (resp. subtract) 1 to the height function at the source to get the 
height function at the target if the edge has a black triangle to the right (resp. left). Also, 
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Figure 1.5 

height functions of a simply connected region R are easy to characterize: height functions 
have the mod 3 values shown in figure 1.5.b, change by at most 2 along any edge in R 
and change by 1 along an edge in the boundary of R. An unexpected consequence of this 
characterization is the following ([STCR]): 

Theorem 1.4: For a simply connected region, let T be the space of its lozenge tilings. 
The set T is a lattice: given two tilings, the maximum (or minimum) of their height 
functions is the height function of some tiling. 

In terms of the pile-of-cubes interpretation, the minimum tiling corresponds to the 
empty box. Thurston ([T]) uses height functions to prove: 

Theorem 1.5: There is a polynomial time algorithm to decide whether a simply con- 
nected region R is tileable by lozenges. 

In a nutshell, start by computing the height function along the boundary of R and 
move inwards assigning to each point the smallest value satisfying the characterization. The 
algorithm either fills the region, obtaining the minimal tiling, or runs into a contradiction, 
indicating non-tileability of R. 

An equivalent criterion, computationally harder but easier to state, is the following: 
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consider positively oriented simple closed curves which, outside the boundary of the region, 
have a black triangle to their left — the numbers of enclosed black and white triangles are 
equal. 

The group theoretical construction easily applies to domino tilings. Label oriented 
edges in the quadriculated plane as in figure 1.1. b, obtaining the Cayley diagram Ti for 
G± = (a, b\[a, b] = e) = Z 2 . Taking into account the shape of the tiles, define G = 
(a,b\[a,b 2 ] = [a 2 ,b] = e): it turns out that G is a non-abelian group and we describe a 
construction of its Cayley diagram V in IR 3 , obtaining a height function for domino tilings. 

Introducing a third generator c = [a, 6], we have G = (a, 6, c\ac = c _1 a, be = c~ 1 b, ab = 
bac) and, by standard arguments, any element of G can be written as c w b y a x for unique 
integers x,y,w. We may thus take Z 3 to be the set of vertices of a Cayley graph for G: 
connect (x, y, w) by an a-edge to (x + 1, y, w) (because c w b y a x ■ a = c w b y a x+1 ), by a 6-edge 
to (x, y + l,w) if x is even or to (x,y + l,w + (— l) y ) if x is odd and by a c-edge to 
(x, y,w + (— l) x+y ). A somewhat nicer picture is obtained by taking 

if x and y are both even, 
if x is odd and y is even, 
if x and y are both odd, 
if x is even and y is odd. 

The vertices of V are then the points (x,y,z) G Z 3 for which 

if x and y are both even, 
if x is odd and y is even, 
if x and y are both odd, 
if x is even and y is odd, 

and, for the original presentation of G, there are edges joining two vertices iff the euclidean 
distance is Notice that this graph projects onto Ti by omitting the z coordinate. 

By construction, any domino tiling lifts to this Cayley graph. In figure 1.4 we show 
a tiling and its related height function. The simple procedure to obtain height functions 
is: colour squares of the region in a checkerboard pattern, assign to the origin and, for 
an edge not covered by a domino, assign to the target vertex the value at the source plus 
(resp. minus) 1 if the square to the left is black (resp. white). 

Again, for simply connected regions, it is easy to give a local characterization of height 
functions, from which the counterparts of theorems 1.4 and 1.5 follow ([STCR] and [T]). 

The pile-of-cubes interpretation for lozenge tilings has its analogue for domino tilings. 
The vertices of the three dimensional solid involved belong to V as above; the solid itself 
can be taken to be 

{(x, y, z) | | a; | < 1, \y\ < 1, |x| + S\y\ — 2\xy\ < z < 4 — 3|x| — \y\ + 2\xy\} , 
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Figure 1.6 



shown in figure 1.7. Notice that, when seen from the top (resp. bottom), the solid covers 
two vertical (resp. horizontal) dominoes. We call this solid a block and there is, indeed, 
a pile-of-blocks interpretation; height functions, however, nicely encapsulate the relevant 
properties of such arrangements and we may therefore avoid visualization altogether. 




Figure 1.7 



The procedures in [CL] and [T] do not apply to regions which are not simply con- 
nected. To study domino tilings, we introduced height sections ([STCR]), a generalization 
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of height functions which applies to any quadriculated surface (e.g., a torus or Klein bot- 
tle). Instead of a formal definition, we give a typical example; vertices in the interior of 
the surface must be surrounded by 4 squares. A similar construction applies for lozenge 
tilings of triangulated surfaces (interior vertices in a triangulated surface are surrounded 
by 6 triangles). 

Consider a cylinder with boundary, as in figure 1.8, tiled by dominoes. The two vertical 
sides are identified as indicated by the labels: two dominoes trespass the cut. Begin by 
assigning the value to the vertex on the left labelled A and follow the local instructions 
for constructing height functions to obtain the value of the height section at all vertices. 
The point A in the cylinder is represented by two vertices in the diagram which received 
different numbers: at the left and 2 at the right. The same ambiguity happens for the 
other points along the cut. How is this to be interpreted? We might take values mod 2 but 
this would amount to throwing away the bundle with the section. Instead, we construct 
a fibre bundle with fibre Z and the cylinder as base space (a 2-manifold with boundary). 
Over any simply connected open subset of the basis, the leaves can be labelled by the 
integers: do this for the complement of the cut — the values of the height section are to be 
understood in terms of these local coordinates. The different values for the points in the 
cut indicate how to glue leaves along the cut: in our example, the n-leaf at the left is glued 
to n + 2-leaf at the right. Each number in the figure labels one element of the fibre over the 
corresponding vertex and the different numbers over distinct representations of the same 
vertex label, by construction, the same element. These elements form a section of the 
restriction of this height bundle to the vertices of the basis: the height section. Summing 
up, the procedure simultaneously constructs the height section and the height bundle. 
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Figure 1.8 



Height sections allow us to give a tileability algorithm similar to Thurston's. It would 
be interesting to detail such an efficient algorithm and generalize the alternate criterion 
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(based on closed curves) to arbitrary quadriculated surfaces. For planar regions, the space 
of tilings is still a lattice. 



2. Counting tilings 

For a quadriculated or triangulated planar region i?, consider a graph A whose vertices 
are the k unit faces (squares or triangles) of R, connected by an edge iff the faces have a 
common side. Since faces in R can be coloured black and white, A is bipartite; we assume 
the number of black faces to equal the number of white faces, a necessary condition for the 
existence of tilings. Let A be the k x k adjacency matrix of A: by labelling black vertices 
first, we may write 




In a language familiar to physicists, tilings of R correspond to dimer coverings of A, i.e., 
partitions of the set of vertices into pairs of adjacent vertices. 

Each dimer covering corresponds to a monomial in the expansion of det B. Suppose 
that B is obtained from B by changing the signs of some entries and that all monomials 
in the expansion of detE have the same sign. Then the number N of tilings of R is 
| detS |. Amazingly, for simply connected triangulated regions, we can take B = B. For 
quadriculated regions, Kasteleyn ([K]) showed that such a B can be obtained from B by 
changing the signs of the l's corresponding to edges on alternate vertical lines. Lieb and 
Loss ([LL]) generalized this idea to prove a formula for the number of dimer partitions of 
an arbitrary bipartite planar graph. The following theorem, a special case of their results, 
unifies the previous claims. 

Theorem 2.1: Let A be a planar bipartite graph with adjacency matrix 




Choose a set of edges in A such that on each minimal loop of size 4n (resp. 4n + 2) there is 
an odd (resp. even) number of chosen edges. Change the sign of the 1 's in B corresponding 
to chosen edges obtaining B. Then the number of dimer partitions of A is | det B |. 

The required choice of edges is always possible. An example is given in figure 2.1, 
where there are 6 minimal loops of sizes 8, 4, 4, 4, 4 and 10 and the chosen edges are thick. 

The sign of the monomial in det B associated to a tiling is not natural, depending 
on the labelling of vertices. Given two tilings, however, equality of their signs is well 
defined and thus tilings naturally split in two classes. This splitting is very different in the 
quadriculated and triangulated cases: 
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Figure 2.1 

Theorem 2.2: Let R be a simply connected region, N the number of tilings of R and B 
as above. If R is quadriculated, 

\detB \ < 1 : 

the two classes are almost balanced. If R is triangulated, 

| det B | = N : 

one of the classes is empty. 

The quadriculated case is due to Deift and Tomei ([DeT]). For general planar regions 
this difference assumes any value. The triangulated case is a corollary of theorem 2.1 but 
we will see another proof in the next section. 

We present closed formulae for N in special cases. Kasteleyn ([K]) obtained a closed 
form for the number of tilings of a rectangle: 

Theorem 2.3: If R is a n x m rectangle, 

m/2 n , j x 1/2 

11 11 V m + 1 n+lj 

k=l£=l v 7 

assuming m even. 

Kasteleyn's proof boils down to computing the determinant of B, as above. He essen- 
tially finds the eigenvalues and eigenvectors of 

(9 B\\ 
\B f J ' 

the eigenvalues are the terms in the product with doubles multiplicities and the eigenvectors 
have the form 

kni £irj 

sin sm . 

m+1 n+1 

The pile-of-cubes interpretation of lozenge tilings reduces their counting for center- 
symmetric hexagonal regions to that of plane partitions. To describe a tiling, count the 
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number of cubes above each square in the xy plane (here, the "up" direction is (0, 0, 1) 
and not (1, 1, 1)). The example in figure 1.3 would be encoded as below. 



5 5 5 5 3 
5 5 110 
4 10 
3 

Following Bender and Knuth ([BK]), a plane partition of n (the total number of cubes in 
the pile) is an array of non-negative integers 

nu n 12 n i3 ••• 
nil n 22 n 23 ■■■ 



for which ^ • • riij = n, > n^+i^j and > n^j+i) (these conditions correspond to 
the pile being gravitationally stable). If = for all i > r (resp. j > c), the partition 
is called r-rowed (resp. c-columned) (this is the size of the region in the x (resp. y) 
direction). If nu < m, we say the parts do not exceed m (the size in the z direction). Let 
7r n be the number of plane r-rowed, c-columned partitions of n with parts not exceeding 
m. MacMahon ([M]) obtained 

2>^ = nn i-^+i-i 

n=0 i=l j—l 

yielding the following formula: 

Theorem 2.4: The number of lozenge tilings of a center-symmetric hexagonal region of 
sides r, c and m is 

r m . . 1 

*=ij=i J 



3. Flips and the space of tilings 

The space T of all tilings of a given region R admits a natural graph structure: two 
tilings are adjacent iff they differ by a flip, the action of lifting two dominoes forming a 
square or three lozenges forming a hexagon and placing them back after a rotation of 90° 
or 60°. Two tilings are joined by a flip iff their height functions are different in one vertex 
only; this difference must be 3 (resp. 4) in the triangulated (resp. quadriculated) case. 
Equivalently, in the pile-of-cubes (or pile-of-blocks) interpretation, a flip corresponds to 
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adding or removing a cube (or block). This suggests a basic property of T when R is 
simply connected: 

Theorem 3.1: If R is a simply connected region, T is connected. 

The inductive step is to prove that if a tiling t\ is larger than a tiling £2, we may 
perform a flip on t\ to obtain ts with t\ > ts > £2- Let $i be the height function for tf, 
remember that all height functions coincide along the boundary. To find out where to flip, 
consider the set of vertices where 9\ — 62 is maximal and, within this set, choose the vertex 
with maximum 6\. Since all height functions coincide in the boundary, this point is a local 
maximum of 6\. By the characterization of height functions, a flip is possible at that point 
(we found a removable cube or block). 

The triangulated version of theorem 2.2 follows from theorem 3.1: for triangulated 
surfaces, flips preserve parity of the monomials in the determinant expansion of B. 

As usual, define distance in the graph T as the minimum number of edges (flips) in a 
path joining two vertices. The following distance formula follows from the same inductive 
step. 

Theorem 3.2: Let R be a simply connected region. The distance between tilings £1 and 
£2, with height functions 9\ and 62, is 

d(*i,*2) = -^X>(t;)-02(v)|, 

v 

where X = 3 or 4 depending whether R is triangulated or quadriculated. 

As usual, we call such minimal paths geodesies. Some geodesies are pretty much the 
same: we may freely change the order of two successive flips ( fi and f'2) if they act on 
disjoint hexagons or squares. We call two such geodesies equivalent and glue a 2-cell (along 
the loop /1/2/r /2~ ) to the graph T to make them homotopic. Similarly, glue fc-cells to 
identify all permutations of k independent flips, obtaining from T the CW-complex T. 

Theorem 3.3: If R is simply connected, T is contractible. 

As a consequence, all geodesies are equivalent and all closed curves null-homotopic. 

We now consider non-simply connected planar regions R. The graph T is now usually 
disconnected: to see this, we consider the following invariant under flips. Take a cut, 
i.e., a simple oriented curve in R consisting of edges with extrema in different boundary 
components of R. If R is triangulated (resp. quadriculated), define the How of a tiling 
across a cut as the number of lozenges (resp. dominoes) crossing the cut, counted with 
a sign determined by the colour of the triangle (resp. square) to the left of the cut. For 
instance, the two tilings in figure 3.1 are not in the same connected component of T because 
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the flows across the indicated cut are 1 and 0. The reader can easily see that homologous 
cuts yield equivalent conditions, the flows differing by a fixed constant from one cut to the 
other. In terms of height sections, it turns out that the height bundle does not depend on 
the tiling and two tilings have the same flow across all cuts iff their height sections coincide 
on the whole boundary. 



Figure 3.1 

Theorem 3.4: Let R be a quadriculated planar region. Two tilings of R belong to the 
same connected component of T iff their Hows coincide across all cuts. For two tilings in the 
same connected component ofT, the distance formula in theorem 3.2 holds. Furthermore, 
connected components of T (constructed as before) are contractible. 

It would be interesting to see what happens for triangulated planar regions. 

A further generalization consists in considering arbitrary triangulated or quadriculated 
surfaces. The cut and flow invariants still make sense if the surface is orientable and 
bicoloured (no longer necessary conditions!). The appropriate generalization to the general 
case involves homology and cohomology with local coefficients, which we shall not discuss. 
In figure 3.2, we show that cut and flow conditions are no longer sufficient, even assuming 
orientability and bicolourability. For obvious reasons, we call these configurations ladders; 
rather surprisingly, their exact position turns out to be the only additional condition for 
domino tilings to be in the same connected component of T. As a further complication, 
if the surface is a (bicolourable) torus, certain connected components of T are no longer 
contractible, but homotopically equivalent to S 1 and therefore contain "closed geodesies" . 
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Theorem 3.5: Let R be an arbitrary quadriculated surface. If the boundary of R 
is non-empty, connected components of T are contractible. In general, they are either 
contractible or homotopically equivalent to circles. 
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